The Rosetta All-Atom Energy Function for Macromolecular Modeling and Design.

Over the past decade, the Rosetta biomolecular modeling suite has informed diverse biological questions and engineering challenges ranging from interpretation of low-resolution structural data to design of nanomaterials, protein therapeutics, and vaccines. Central to Rosetta's success is the energy function: a model parametrized from small-molecule and X-ray crystal structure data used to approximate the energy associated with each biomolecule conformation. This paper describes the mathematical models and physical concepts that underlie the latest Rosetta energy function, called the Rosetta Energy Function 2015 (REF15). Applying these concepts, we explain how to use Rosetta energies to identify and analyze the features of biomolecular models. Finally, we discuss the latest advances in the energy function that extend its capabilities from soluble proteins to also include membrane proteins, peptides containing noncanonical amino acids, small molecules, carbohydrates, nucleic acids, and other macromolecules.

Biological function derived from predicted structures in CASP11.

In CASP11, the organizers sought to bring the biological inferences from predicted structures to the fore. To accomplish this, we assessed the models for their ability to perform quantifiable tasks related to biological function. First, for 10 targets that were probable homodimers, we measured the accuracy of docking the models into homodimers as a function of GDT-TS of the monomers, which produced characteristic L-shaped plots. At low GDT-TS, none of the models could be docked correctly as homodimers. Above GDT-TS of ∼60%, some models formed correct homodimers in one of the largest docked clusters, while many other models at the same values of GDT-TS did not. Docking was more successful when many of the templates shared the same homodimer. Second, we docked a ligand from an experimental structure into each of the models of one of the targets. Docking to the models with two different programs produced poor ligand RMSDs with the experimental structure. Measures that evaluated similarity of contacts were reasonable for some of the models, although there was not a significant correlation with model accuracy. Finally, we assessed whether models would be useful in predicting the phenotypes of missense mutations in three human targets by comparing features calculated from the models with those calculated from the experimental structures. The models were successful in reproducing accessible surface areas but there was little correlation of model accuracy with calculation of FoldX evaluation of the change in free energy between the wild-type and the mutant. Proteins 2016; 84(Suppl 1):370-391. © 2016 Wiley Periodicals, Inc.

Stability engineering of anti-EGFR scFv antibodies by rational design of a lambda-to-kappa swap of the VL framework using a structure-guided approach.

Phage-display technology facilitates rapid selection of antigen-specific single-chain variable fragment (scFv) antibodies from large recombinant libraries. ScFv antibodies, composed of a VH and VL domain, are readily engineered into multimeric formats for the development of diagnostics and targeted therapies. However, the recombinant nature of the selection strategy can result in VH and VL domains with sub-optimal biophysical properties, such as reduced thermodynamic stability and enhanced aggregation propensity, which lead to poor production and limited application. We found that the C10 anti-epidermal growth factor receptor (EGFR) scFv, and its affinity mutant, P2224, exhibit weak production from E. coli. Interestingly, these scFv contain a fusion of lambda3 and lambda1 V-region (LV3 and LV1) genes, most likely the result of a PCR aberration during library construction. To enhance the biophysical properties of these scFvs, we utilized a structure-based approach to replace and redesign the pre-existing framework of the VL domain to one that best pairs with the existing VH. We describe a method to exchange lambda sequences with a more stable kappa3 framework (KV3) within the VL domain that incorporates the original lambda DE-loop. The resulting scFvs, C10KV3_LV1DE and P2224KV3_LV1DE, are more thermodynamically stable and easier to produce from bacterial culture. Additionally, C10KV3_LV1DE and P2224KV3_LV1DE retain binding affinity to EGFR, suggesting that such a dramatic framework swap does not significantly affect scFv binding. We provide here a novel strategy for redesigning the light chain of problematic scFvs to enhance their stability and therapeutic applicability.

BioAssemblyModeler (BAM): user-friendly homology modeling of protein homo- and heterooligomers.

UNLABELLED: Many if not most proteins function in oligomeric assemblies of one or more protein sequences. The Protein Data Bank provides coordinates for biological assemblies for each entry, at least 60% of which are dimers or larger assemblies. BioAssemblyModeler (BAM) is a graphical user interface to the basic steps in homology modeling of protein homooligomers and heterooligomers from the biological assemblies provided in the PDB. BAM takes as input up to six different protein sequences and begins by assigning Pfam domains to the target sequences. The program utilizes a complete assignment of Pfam domains to sequences in the PDB, PDBfam (http://dunbrack2.fccc.edu/protcid/pdbfam), to obtain templates that contain any or all of the domains assigned to the target sequence(s). The contents of the biological assemblies of potential templates are provided, and alignments of the target sequences to the templates are produced with a profile-profile alignment algorithm. BAM provides for visual examination and mouse-editing of the alignments supported by target and template secondary structure information and a 3D viewer of the template biological assembly. Side-chain coordinates for a model of the biological assembly are built with the program SCWRL4. A built-in protocol navigation system guides the user through all stages of homology modeling from input sequences to a three-dimensional model of the target complex. AVAILABILITY: http://dunbrack.fccc.edu/BAM.

Minimal ensembles of side chain conformers for modeling protein-protein interactions.

The goal of this article is to reduce the complexity of the side chain search within docking problems. We apply six methods of generating side chain conformers to unbound protein structures and determine their ability of obtaining the bound conformation in small ensembles of conformers. Methods are evaluated in terms of the positions of side chain end groups. Results for 68 protein complexes yield two important observations. First, the end-group positions change less than 1 Å on association for over 60% of interface side chains. Thus, the unbound protein structure carries substantial information about the side chains in the bound state, and the inclusion of the unbound conformation into the ensemble of conformers is very beneficial. Second, considering each surface side chain separately in its protein environment, small ensembles of low-energy states include the bound conformation for a large fraction of side chains. In particular, the ensemble consisting of the unbound conformation and the two highest probability predicted conformers includes the bound conformer with an accuracy of 1 Å for 78% of interface side chains. As more than 60% of the interface side chains have only one conformer and many others only a few, these ensembles of low-energy states substantially reduce the complexity of side chain search in docking problems. This approach was already used for finding pockets in protein-protein interfaces that can bind small molecules to potentially disrupt protein-protein interactions. Side-chain search with the reduced search space will also be incorporated into protein docking algorithms.

Nonplanar peptide bonds in proteins are common and conserved but not biased toward active sites.

The planarity of peptide bonds is an assumption that underlies decades of theoretical modeling of proteins. Peptide bonds strongly deviating from planarity are considered very rare features of protein structure that occur for functional reasons. Here, empirical analyses of atomic-resolution protein structures reveal that trans peptide groups can vary by more than 25° from planarity and that the true extent of nonplanarity is underestimated even in 1.2 Å resolution structures. Analyses as a function of the ϕ,ψ-backbone dihedral angles show that the expected value deviates by ± 8° from planar as a systematic function of conformation, but that the large majority of variation in planarity depends on tertiary effects. Furthermore, we show that those peptide bonds in proteins that are most nonplanar, deviating by over 20° from planarity, are not strongly associated with active sites. Instead, highly nonplanar peptides are simply integral components of protein structure related to local and tertiary structural features that tend to be conserved among homologs. To account for the systematic ϕ,ψ-dependent component of nonplanarity, we present a conformation-dependent library that can be used in crystallographic refinement and predictive protein modeling.

A smoothed backbone-dependent rotamer library for proteins derived from adaptive kernel density estimates and regressions.

Rotamer libraries are used in protein structure determination, prediction, and design. The backbone-dependent rotamer library consists of rotamer frequencies, mean dihedral angles, and variances as a function of the backbone dihedral angles. Structure prediction and design methods that employ backbone flexibility would strongly benefit from smoothly varying probabilities and angles. A new version of the backbone-dependent rotamer library has been developed using adaptive kernel density estimates for the rotamer frequencies and adaptive kernel regression for the mean dihedral angles and variances. This formulation allows for evaluation of the rotamer probabilities, mean angles, and variances as a smooth and continuous function of phi and psi. Continuous probability density estimates for the nonrotameric degrees of freedom of amides, carboxylates, and aromatic side chains have been modeled as a function of the backbone dihedrals and rotamers of the remaining degrees of freedom. New backbone-dependent rotamer libraries at varying levels of smoothing are available from http://dunbrack.fccc.edu.

Development of a rotamer library for use in beta-peptide foldamer computational design.

Foldamers present a particularly difficult challenge for accurate computational design compared to the case for conventional peptide and protein design due to the lack of a large body of structural data to allow parametrization of rotamer libraries and energies. We therefore explored the use of molecular mechanics for constructing rotamer libraries for non-natural foldamer backbones. We first evaluated the accuracy of molecular mechanics (MM) for the prediction of rotamer probability distributions in the crystal structures of proteins is explored. The van der Waals radius, dielectric constant and effective Boltzmann temperature were systematically varied to maximize agreement with experimental data. Boltzmann-weighted probabilities from these molecular mechanics energies compare well with database-derived probabilities for both an idealized alpha-helix (R = 0.95) as well as beta-strand conformations (R = 0.92). Based on these parameters, de novo rotamer probabilities for secondary structures of peptides built from beta-amino acids were determined. To limit computational complexity, it is useful to establish a residue-specific criterion for excluding rare, high-energy rotamers from the library. This is accomplished by including only those rotamers with probability greater than a given threshold (e.g., 10%) of the random value, defined as 1/n where n is the number of potential rotamers for each residue type.

Conformation dependence of backbone geometry in proteins.

Protein structure determination and predictive modeling have long been guided by the paradigm that the peptide backbone has a single, context-independent ideal geometry. Both quantum-mechanics calculations and empirical analyses have shown this is an incorrect simplification in that backbone covalent geometry actually varies systematically as a function of the phi and Psi backbone dihedral angles. Here, we use a nonredundant set of ultrahigh-resolution protein structures to define these conformation-dependent variations. The trends have a rational, structural basis that can be explained by avoidance of atomic clashes or optimization of favorable electrostatic interactions. To facilitate adoption of this paradigm, we have created a conformation-dependent library of covalent bond lengths and bond angles and shown that it has improved accuracy over existing methods without any additional variables to optimize. Protein structures derived from crystallographic refinement and predictive modeling both stand to benefit from incorporation of the paradigm.

Improved prediction of protein side-chain conformations with SCWRL4.

Determination of side-chain conformations is an important step in protein structure prediction and protein design. Many such methods have been presented, although only a small number are in widespread use. SCWRL is one such method, and the SCWRL3 program (2003) has remained popular because of its speed, accuracy, and ease-of-use for the purpose of homology modeling. However, higher accuracy at comparable speed is desirable. This has been achieved in a new program SCWRL4 through: (1) a new backbone-dependent rotamer library based on kernel density estimates; (2) averaging over samples of conformations about the positions in the rotamer library; (3) a fast anisotropic hydrogen bonding function; (4) a short-range, soft van der Waals atom-atom interaction potential; (5) fast collision detection using k-discrete oriented polytopes; (6) a tree decomposition algorithm to solve the combinatorial problem; and (7) optimization of all parameters by determining the interaction graph within the crystal environment using symmetry operators of the crystallographic space group. Accuracies as a function of electron density of the side chains demonstrate that side chains with higher electron density are easier to predict than those with low-electron density and presumed conformational disorder. For a testing set of 379 proteins, 86% of chi(1) angles and 75% of chi(1+2) angles are predicted correctly within 40 degrees of the X-ray positions. Among side chains with higher electron density (25-100th percentile), these numbers rise to 89 and 80%. The new program maintains its simple command-line interface, designed for homology modeling, and is now available as a dynamic-linked library for incorporation into other software programs.

Vancomycin forms ligand-mediated supramolecular complexes.

The emergence of resistance to vancomycin and related glycopeptide antibiotics is spurring efforts to develop new antimicrobial therapeutics. High-resolution structural information about antibiotic-ligand recognition should prove valuable in the rational design of improved drugs. We have determined the X-ray crystal structure of the complex of vancomycin with N-acetyl-D-Ala-D-Ala, a mimic of the natural muramyl peptide target, and refined this structure at a resolution of 1.3 A to R and R(free) values of 0.172 and 0.195, respectively. The crystal asymmetric unit contains three back-back vancomycin dimers; two of these dimers participate in ligand-mediated face-face interactions that produce an infinite chain of molecules running throughout the crystal. The third dimer packs against the side of a face-face interface in a tight "side-side" interaction that involves both polar contacts and burial of hydrophobic surface. The trimer of dimers found in the asymmetric unit is essentially identical to complexes seen in three other crystal structures of glycopeptide antibiotics complexed with peptide ligands. These four structures are derived from crystals belonging to different space groups, suggesting that the trimer of dimers may not be simply a crystal packing artifact and prompting us to ask if ligand-mediated oligomerization could be observed in solution. Using size-exclusion chromatography, dynamic light scattering, and small-angle X-ray scattering, we demonstrate that vancomycin forms discrete supramolecular complexes in the presence of tripeptide ligands. Size estimates for these complexes are consistent with assemblies containing four to six vancomycin monomers.

BioDownloader: bioinformatics downloads and updates in a few clicks.

UNLABELLED: There are many ftp or http servers storing data required for biological research. While some download applications are available, there is no user-friendly download application with a graphical interface specifically designed and adapted to meet the requirements of bioinformatics. BioDownloader is a program for downloading and updating files from ftp and http servers. It is optimized to work robustly with large numbers of files. It allows the selective retrieval of only the required files (batch downloads, multiple file masks, ls-lR file parsing, recursive search, recent updates, etc.). BioDownloader has a built-in repository containing the settings for common bioinformatics file-synchronization needs, including the Protein Data Bank (PDB) and National Center for Biotechnology Information (NCBI) databases. It can post-process downloaded files, including archive extraction and file conversions. AVAILABILITY: The program can be installed from http://dunbrack.fccc.edu/BioDownloader. The software is freely available for both non-commercial and commercial users under the BSD license.

Statistical and conformational analysis of the electron density of protein side chains.

Protein side chains make most of the specific contacts between proteins and other molecules, and their conformational properties have been studied for many years. These properties have been analyzed primarily in the form of rotamer libraries, which cluster the observed conformations into groups and provide frequencies and average dihedral angles for these groups. In recent years, these libraries have improved with higher resolution structures and using various criteria such as high thermal factors to eliminate side chains that may be misplaced within the crystallographic model coordinates. Many of these side chains have highly non-rotameric dihedral angles. The origin of side chains with high B-factors and/or with non-rotameric dihedral angles is of interest in the determination of protein structures and in assessing the prediction of side chain conformations. In this paper, using a statistical analysis of the electron density of a large set of proteins, it is shown that: (1) most non-rotameric side chains have low electron density compared to rotameric side chains; (2) up to 15% of chi1 non-rotameric side chains in PDB models can clearly be fit to density at a single rotameric conformation and in some cases multiple rotameric conformations; (3) a further 47% of non-rotameric side chains have highly dispersed electron density, indicating potentially interconverting rotameric conformations; (4) the entropy of these side chains is close to that of side chains annotated as having more than one chi(1) rotamer in the crystallographic model; (5) many rotameric side chains with high entropy clearly show multiple conformations that are not annotated in the crystallographic model. These results indicate that modeling of side chains alternating between rotamers in the electron density is important and needs further improvement, both in structure determination and in structure prediction.